Two-eigenfunction correlation in a multifractal metal and insulator 
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We consider the correlation of two single-particle probability densities |4'_B(r)p at coinciding 
points r as a function of the energy separation ui = \E—E'\ for disordered tight-binding lattice models 
(the Anderson models) and certain random matrix ensembles. We focus on the parameter range close 
but not exactly at the Anderson localization transition. We show that even away from the critical 
point the eigenfunction statistics exhibit the remnant of multifractality characteristic of the critical 
states. This leads to an enhancement of eigenfunction correlations and a corresponding enhancement 
of matrix elements of the local electron interaction at small energy separations. This enhancement 
is accompanied by a depression of correlations at large energy separations, both phenomena being 
a consequence of the stratification of space into densely packed but mutually avoiding resonance 
clusters. We also demonstrate that the correlation function of localized states in a d-dimensional 
insulator is logarithmically enhanced at small energy separations provided that d > 1. A simple and 
general physical picture of all these phenomena is presented. 

Finally by a combination of numerical results on the Anderson model and analytical and numerical 
results for the relevant random matrix theories we identified the Gaussian random matrix ensembles 
that describe the multifractal features both in the metal and in the insulator phases. 

PACS numbers: 72.15.Rn, 72.70.+m, 72.20.Ht, 73.23.-b 
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I. INTRODUCTION 

Eigenfunction and spectral statistics in quantum sys- 
tems with quenched disorder were a subject of intense 
studj^ in the context of mesoscopic fluctuations of con- 
ductance and density of states (DoS), in particular in 
quantum dots^ . For this application the most relevant is 
the regime of weak deviation^i^ from the Wigner-Dyson 
statistics given by the conventional random matrix the- 
ory (RMT)^. Disordered multi-channel quantum wires 
is the most important example of systems where single- 
particle eigenstates are all localized. Here the statistics 
of eigenstates require a non-perturbative treatment us- 
ing the formalism of nonlinear sigma-model^ or banded 
random matrices^. A special class are systems with the 
critical, multifractal (MF) eigenstate statisticsii^i^. Two- 
dimensional disordered metals fall in this class-- provided 
that effects of localization are suppressed by magnetic 
field. Otherwise, one can speak only on weak multi- 
fractality which turns to localization before being fully 
developed. The true physical realizations of the crit- 
ical, MF eigenstate statistics are systems at the crit- 
ical point of the Anderson localization transitio n^^i^^ 
and the integer Quantum Hall systems at the center of 
Landau bandit. Importantly, the class of systems with 
MF eigenstate statistics also allows for a random ma- 
trix representation^, in particular using the power-law 
banded random matrices (PLBRM)^i. 

Another field of intense research is the interplay be- 
tween disorder and electron interaction with the seminal 
results on quantum correction to the tunnel DoS and 
conductivitj^ of disordered two-dimensional metals and 



the correction to superconducting transition temperature 
due to a simultaneous effect of disorder and the Coulomb 
interactioni^. In all those works disorder and interaction 
are taken into account essentially perturbatively along 
the lines given inii. Recently there was an attempti^ to 
consider the problem of superconductivity near the An- 
derson transition in which disorder has been treated non- 
perturbatively by postulating the MF statistics of one- 
particle states ^i(r) that enter the matrix element of a 
phenomenological electron attraction: 



Ji] = 9 



J dr^,(r)2*j (r)2. 



(1) 



In particular, the simplest quantity of interest is the dis- 
order average matrix element (Jy ) at a given energy sep- 
aration to between one-particle energies Ei and Ej. For 
real eigenf unctions (orthogonal symmetry class) it is pro- 
portional to the correlation function C(w) = K{uj)/R{lu), 
where 



dr 



(\^,iv)\^\^,{r)\HiE,-E,^u)), (2) 



2J 



and R{uj) = J2ij{^i-^i ^ ~ ^)) the spectral corre- 
lation function which is close to 1 for uj much greater 
then the mean level spacing. The correlation function 
C(w) ~ K(lu) is the main subject of the present paper. 

The correlation function defined by Eq.® is a mea- 
sure of overlap of two different eigenfunctions. For truly 
extended normalized states (e.g. in a quantum dot) 
|*i(r)p = 1/V and thus C{uj)V = 1, where V is a sys- 
tem volume. Remarkably, C{uj)V — 1 is also valid for 
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classical examples of localized states, e.g. in a quantum 
disordered wire. In this case two states are typically not 
overlapping but with a small probability of ^''/V (where 
^ is the localization radius) they are localized in the same 
place and then the integral in Eq.(l2]) is of the order of 
the inverse localization volume 1/^'*. 

There are cases, however, when an eigenfunction ^'i(r) 
does not occupy all the available volume or all the local- 
ization volume and the typical amplitude |5'i(r)p is not 
just the inverse volume (for extended states) or the in- 
verse localization volume (for localized states). In this 
case a non-trivial behavior of the correlation function 
C{uj) is expected. Such situation is realized near the 
critical point of the Anderson localization transition. In 
the vicinity of this point in the region of extended states 
{multifractal metal) or in the region of localized states 
{multifractal insulator) the system retains the character- 
istic features of the critical multifractal statistics of eigen- 
states which makes it qualitatively different from both a 
normal metal or a normal Anderson insulator. 

In this paper we will identify and quantify such char- 
acteristic features in the correlation function C{uj) and 
give their interpretation in terms of the typical behav- 
ior of single-particle states. To attain this goal we will 
combine new analytical results for the PLBRM with nu- 
merics on the PLBRM and the Anderson model. We 
specially focus on the dependence of C{u!) on the energy 
difference to in the cross-over region in the vicinity but 
not exactly at the Anderson transition point which has 
not been studied so far. 

The paper is organized as follows. In section II we give 
a brief introduction into the subject of multifractality of 
critical eigenstates focusing on the main effect of multi- 
fractality which is the critical enhancement of eigenfunc- 
tion correlations. In section III we give a cartoon of the 
off-critical states in a multifractal metal and a multifrac- 
tal insulator and introduce the random matrix theories 
which may describe them. In Sec. IV we present the re- 
sults of an analytical theory of eigenfunction correlations 
for a class of almost diagonal Gaussian random matri- 
ces which all the RMT's suggested to describe strong 
multifractality fall into. In Sec.V we consider the two- 
eigenfunction correlation function exactly at the critical 
point of the localization transition in the 3D Anderson 
model and for the critical random matrix ensemble in 
the limit of strong multifractality. We show that the dy- 
namical scaling relationship suggested by Chalker is not 
violated even in the limit when the fractal dimensions 
are very small. In Section VI we describe the new phe- 
nomenon of eigenfunction mutual avoiding and present 
a qualitative picture that simultaneously explains the 
enhancement of eigenfunction correlations at small en- 
ergy separations and the eigenfunction mutual avoiding 
at large energy separations. In section VII we consider 
the properties of eigenfunction correlations in a multi- 
fractal insulator. In particular, we describe the new phe- 
nomenon of logarithmic enhancement of eigenfunction 
correlations at small energy separations in the 2D and 



3D Anderson insulators and show the absence of such en- 
hancement in the quasi-lD case. We also suggest a Trun- 
cated Critical RM ensemble that describes all the princi- 
pal features of eigenfunction correlations in the 3D multi- 
fractal insulator. Section VIII is devoted to the random 
matrix description of the multifractal metal. We show 
that the sub-critical PLBRM suggested in Ref.^'* gives a 
reasonable agreement with the 3D Anderson model. By 
analytical treatment of this RM model we found the re- 
gion of parameters where the eigenfunction correlations 
become effectively short-range in the energy space which 
may poit out on the existence of a new metal phase above 
some critical dimensionlity in the multi-dimensional An- 
derson model. In the Conclusion we list all the principal 
results of this paper. 



II. MULTIFRACTALITY OF CRITICAL 
EIGENFUNCTIONS 

The "standard" model (the Anderson model) for the 
Anderson localization transition in d > 2 dimensions is 
the tight-binding model with the hopping constant V — 1 
and random on-site energies Ei characterized by the dis- 
tribution function V{si) which is frequently chosen con- 
stant V{ei) = 1/W in the interval [~W/2, W/2] and zero 
otherwise. There is a vast literature (see e.g^ and ref- 
erences therein) on numerical investigation of the Ander- 
son localization transition in this model on a 3D lattice. 
Recently also higher dimensions d > 3 become accessi- 
ble to modern computers^S. While the earlier studies of 
this model were focused on the critical behavior of the 
localization/correlation length ^ near the critical disor- 
der Wc, the recent works were mostly related with the 
statistics of critical eigenfunctions. The multifractality 
of critical eigenfunctions predicted inii almost immedi- 
ately after emergence of scaling theory of localization has 
been confirmed and quantified in detail. 

The results obtained for the Anderson model exactly 
at the critical point seem to be very well described^ ^ by 
the critical PLBRM model^ii^. This model is defined as 
an ensemble of random Hermitean matrices which en- 
tries Hij fluctuate independently around zero with the 
variance: 

where (3 = 1,2,4 for the Dyson orthogonal, unitary, 
and symlectic symmetry classes^ and b is the parameter 
that controls the multifractality exponents. This model 
has been studied and its comparison with the Anderson 
model in d dimensions has been done predominantly for 
the statistical moments P„ of a single eigenstate at a 
given energy E: 

P,{E) = ^ ^(|*„(r)p' 5{E E„)). (4) 
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The best known example is the inverse participation ratio 
(IPR) given by the second moment P2 . The multifractal 
statistics of a single eigenstate is characterized by the 
moment Pq that scales with the system volume V or the 
total number of sites N as: 



(5) 



where dq < d is the fractal dimension corresponding to 
q-th moment. The existence of the scaling law Eq.([5|) and 
the dependence of the exponent dq on q are the principle 
features of eigenfunction multifractality. The fractal ex- 
ponents dq depend also on the symmetry class f3 and the 
space dimensionality d. For the critical PLBRM Eq.Q 
the dependence on d is modeled by the dependence of dq 
on the parameter b. 

The critical scaling Eq.(l5]) with respect to the system 
size N has its dynamical counterpart when instead of one 
single eigenfunction one considers the correlation func- 
tion Eq. ([2|) of two eigenfunctions at an energy separation 
u! = \E — E'\ between them. This scaling has been sug- 
gested by Chalker— many years ago: 



S <iLj < Eo, 



where 



1 ^2 



(6) 



(7) 



S is the mean level spacing and Eq is the upper cut- 
off of multifractality. Numerics on the integer Quantum 
Hall systems and in the critical point of the 3D Anderson 
model was consisten t^^i^^ with this scaling. 

An important feature of Eq.® is that the exponent 
1 — ^2/1^ in the w-dependence is smaller than 1. Even in 
the limit of infinitely small correlation dimension ^2 the 
correlation function decays slowly as 1/lo. This implies 
that the sparse critical states separated by large energy 
distance are still well overlaping^^, in contrast to strongly 
localized states which typically do not overlap even for 
nearest neighbors in the energy space. The reason for 
such a behavior and the physical meaning of the energy 
scale Eq will be discussed in Section VI. 

As the correlation function is equal to C(w) ~ 
both for the truly extended and the ideal localized states, 
Eq. ^ implies the critical enhancement of eigenfunction 
correlations for lo < Eq. This enhancement is cru- 
cially important for electron interaction near the An- 
derson localization transition, in particular for the su- 
perconducting transition temperature^®. To illustrate 
this point we present in Fig.l the result of numeri- 
cal diagonalization of the critical PLBRM, the classi- 
cal Wigner-Dyson RM (jiJijp) = const with extended 
eigenstates, and the ensemble of conventional banded 
random matrices^ with exponentially decreasing entries 
(l-^iiP) exp{—\i — j\/B) which describes strongly lo- 
calized eigenstates in quasi-lD disordered systems. The 
critical enhancement of eigenfunction correlations is evi- 
dent from this plot. 
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FIG. 1: Critical enhancement of eigenfunction correlation. 
Results of exact diagonalization of the critical PLBRM at 
b = 0.1, the banded random matrices with B=5, and Wigner- 
Dyson RM are shown in red, blue and green, respectively, 
and are represented by squares (N=200), circles (N=1000) 
and stars (N=2000). 



The physical origin of the enhancement is two-fold: 
(i) a critical eigenfunction "occupies" only part of the 
available space which by normalization — 1 

enhances its amplitude, and (ii) the supports (the man- 
ifold of {r} where |\I'(r)p is essentially non-zero) of dif- 
ferent critical eigenfunctions are strongly overlapping. It 
is important that both conditions are fulfilled simultane- 
ously. For instance the condition (i) is fulfilled for local- 
ized states even better than for the critical ones but the 
lack of the condition (ii) levels off the gain in the cor- 
relation function C{uj). On the contrary, in a metal the 
condition (ii) is trivially fulfilled, but the eigenfunction 
amplitude is small. 



III. OFF-CRITICAL STATES AND THEIR 
RANDOM MATRIX REPRESENTATIONS 



Gaussian random matrix models proved to be an effi- 
cient and universal theoretical tool for describing com- 
plex systems. The success was partially due to the 
available analytical solutions^ii and partially due to effi- 
cient algorithms of numerical diagonalization of matrices. 
Therefore it is highly desirable to have random matrix 
models that describe not only the critical MF eigenstates 
but also localized and extended eigenstates in the vicin- 
ity of the Anderson transition. The criterion to select 
such models is a qualitative and (when possible) a quan- 
titative agreement with the results on the 3D Anderson 
model. 
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FIG. 2: 2D cartoon of a) conventional localized state; b) lo- 
calized state in a multifractal insulator; c) extended state in a 
multifractal metal. The darker regions correspond to higher 
eigenfunction amplitude. The localization/correlation radius 
^ is shown in each case. 



As will be demonstrated below, the correlation func- 
tion C{uo) in the 3D Anderson model contains the crit- 
ical power-law behavior Eq.Q well beyond the Ander- 
son transition point. As a matter of fact the correlation 
function C{lj) is indistinguishable from the critical one 
until the dynamic length = l/{pujY/'^ (p is the mean 
DoS) exceeds the localization/correlation length ^. For 
Luj > or u smaller than the level spacing in the local- 
ized volume ^ l/{p£f')^ the correlation function loses 
its critical features and shows typical features of a metal 
or an isulator. 

This allows us to suggest the following cartoon of typi- 
cal eigenf unctions in the vicinity of the localization tran- 
sition shown in Fig. 2. Namely, a typical localized state 
in a "multifractal insulator" can be viewed as a "piece 
of multifractal" of the size of the localization radius ^ 
(Fig. 2b.), in contrast to a conventional localized state 
where all the localization volume is more of less homo- 
geneously "filled" (Fig. 2a.). In the same way, typical 
extended states on the metallic side of the localization 
transition ( "multifractal metal") should look like a mo- 
saic made of such "pieces of multifractal" (Fig. 2c.). 

Based on the persistence of the critical behavior be- 
yond the critical region it is natural to assume that the 
random matrix model for the extended states near the 
critical point and the localized states on the other side of 
the transition should bear features of the critical RMT. 

Let us start by constructing a random matrix model 
for the multifractal insulator. Given that the quasi- 
ID insulator is well described by the banded ran- 
dom matrices^, with exponentially decaying variance 
{\Hij\^) = exp{— |i — i\/B} we suggest the following hy- 



brid RM model as a model for the multifractal insulator: 

^^1 2(1.-,%-) exp{-(|»-j|/m, 

(8) 

As compared with the critical PLBRM model Eq.©, the 
model Eq.Q contains an additional parameter B which 
sets in a finite localization radius £,{B). It also contains 
an exponent ry which depends on the space dimensionality 
d of the disordered lattice model we would like to model 
by the RMT. In Sec. VIII we give both analytical and 
numerical arguments in favor of the choice 

7J = 1/d. 

Another candidate has been suggested in Ref J^: 

In this case the localization radius ^ is controlled by the 
variable exponent a of the power-law. For a multifractal 
insulator a > 1. 

The possible RM models for multifractal metal are also 
constructed as deformations of the critical PLBRM. The 
model Eq. ([9]) for a < 1 is believedi^ to describe the mul- 
tifractal metal. One can also think that the Gaussian 
RMT 

which is a hybrid of the critical PLBRM and the WD 
RMT, is also suitable for this purpose. Below we will 
study all those RM models in detail and compare the 
corresponding results for the correlation function C(w) 
with the results obtained by numerical diagonalization 
of the d-dimensional Anderson model. 



IV. ALMOST DIAGONAL GAUSSIAN RMT: 
ANALYTICAL RESULTS FOR C{lu) 

The characteristic properties of multifractal statistics 
of critical and off-critical states are best seen when the 
multifractality is strong. This is the case where the pa- 
rameter & < 1 in Eqs.jSl m m [To]) is small. On the 
other hand, this is exactly the limit where the typical 
off-diagonal elements of Hij are small compared to di- 
agonal ones. Such matrices (referred to as almost diago- 
nal raxidoiji matrices (ADRM)) may possess a non-trivial 
statistics of eigenfunctions which justifies their special 
study^i^. The idea of analytical treatment of ADRM, 
first suggested in Refi^^ and used in^l to compute the 
correlation dimension d2 for the critical PLBRM model 
Eq.Q, is similar to the virial expansion in dilute gases. 
However, instead of taking into account two-, three- and 
multiple-particle collisions, one considers progressively 
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increasing number of interacting resonance sites coupled 
by a small off-diagonal matrix element H^j . Recently the 
virial coefficients for the Gaussian ADRM with an ar- 
bitrary (but small) variance (t^(|i — j|) = 2(|i7yp) were 
expressed through the supersymmetric field theory^^ and 
the correlation function C{lu) has been explicitly calcu- 
lated in the two-state approximation for the unitary sym- 
metry class (3 = 2. The result is the following: 



kjuj) 



where 



N 



n=l 



2u)+ - 



Erfi (uj) ~ 1 
2 ^ ' 



and 



''H-^E'^^"^'Erfi(<:^). (13) 

n=l 



In Eqs. (|12ll3p we denote 

Ld LjJ 



(14) 



and Erfi(z) = ^ di. 

The result given by Eas. p2ll3p is valid in the limit 
when: 



N 



(15) 



For the RMT defined by Eq.® and Eq.® with a> I 
which are suggested to describe the multifractal insula- 
tor, the sum over n in Eq. ([15]) converges. Then the va- 
lidity of Eqs. p2ll3p is independent of the matrix size TV 
in the limit N oo and is controlled only by a small 
parameter & <C 1. On the contrary, for the models of the 
multifractal metal described by Eq. p^ and Eq.® with 
a < 1 the sum in Eq.ljTS]) diverges at large N. Then 
Eqs. (|12ll3p arc only valid for iV < ^ where we define the 
correlation radius ^ as follows: 



(16) 



n=l 



We will show below that a good qualitative description 
of the metal phase in the limit N ^ oo can still be ob- 
tained from the above theory if one substitutes ^ for N 
in Eqs.dlini). 

Eqs. p2l) - p6l) will be used throughout the paper to 
analize different random matrix ensembles suggested as 
possible models for critical eigenstates and the off-critical 
states in a multifractal metal and insulator. 



TWO EIGENFUNCTION CORRELATIONS 
AT CRITICALITY 



It is not a priori clear that the critical power-law behav- 
ior Eq.® and the dynamical scaling relationship Eq.Q 
hold true for all systems where Eq.([5]) is valid. In par- 
ticular it is interesting to study the correlation function 
Eq.(I2]) in the limit of strong multifractality when ^2^0. 
Below we will derive an analytical formula for the criti- 

(11) cal PLBRM in the limit b ^ which corresponds^" to 
d2 and confirm the scaling law Eq.([7]) by numerical 
diagonalization of PLBRM with very small b. 

One can easily see from Eqs. (IT^ - p^ in which we plug 
in Eq.® that in the interval < \E - E'\ < 6 the 

(12) correlation function C{oj) given by these equations has 
an asymptotic power law behavior Eq.® with /.t = 1: 



\lu\ 



'7r\3/2 

£^0 = b. (17) 



Applying Eq.HS]) to the critical PLBRM Eq.jHl) gives 
the criterion of validity InA^ <C 1/b ^ l/d2 in which 
case one cannot distinguish between N-^"'^'^ and N, or, 
correspondingly, between uj^^^'^^ and l/oj. Thus the an- 
alytical formulae Eqs. (|12l3p is consistent with the scaling 
relationship Eq.(l7]), given that d2{b) ^ as ^ 0. 

A comparison of the analytical results and the results 
of numerical diagonalization of the critical PLBRM with 
(3 = 2 and b = 0.06 is shown in Fig. 3. The coincidence is 
very good for large energy separations. The deviation at 
small energy separations is due to the difference in the 
values of /x. The exponent = 1 for the analytical curve 
and n = 0.86 for the numerical curve which is very close 
to the prediction of Eq.([7]) /i = 1 — ~ 0.865, where d2 
is found from the numerical data for P2{N) and Eq.®. 

The scaling relationship Eq.([7]) is further checked in 
Fig. 4 where the numerical data for fi and 1 — (i2 is plot- 
ted as a function of b. The fulfillment of this relationship 
down to b as small as 0.005 and an agreement with the 
theoretical prediction of Refill is spectacular. Thus 
from the combination of analytical and numerical results 
we conclude that the Chalker's scaling Eqs.®,© is valid 
for arbitrary small b and thus for arbitrary strong multi- 
fractality. 

Finally we demonstrate how well the critical PLBRM 
Eq. ^ describes the two-eigenfunction correlations in the 
3D Anderson model at the mobility edge. To this end we 
modify the distribution of the on-site energies in the An- 
derson model from the standard rectangular box distri- 
bution to the triangular distribution where the mobility 
edge corresponds to Ec = ±3.5. The correlation func- 
tion C{u!) with E, E' near the mobility edge is shown in 
Fig. 5. It coincides almost exactly with the correspond- 
ing curve resulting from numerical diagonalization of the 
critical PLBRM ensemble with only one fitting parame- 
ter b = 0.42. 
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FIG. 3: Two-eigenfunction correlation function C{uj) for 
the critical PLBRM with /3 = 2, 6 = 0.06 and N = 
200(square),1000(circle) and 2000 (star). The analytical 
curve at A*' = 2000 given by Eqs. (|llj|T5)) is shown by a solid 
line. 
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FIG. 4: The scaling relationship between /i and d2 for the 
/3 = 2 critical PLBRM. The solid line is the prediction based 
on Ref.^ d2 = ^b. 
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VI. EIGEFUNCTION MUTUAL AVOIDING 
AND STRATIFICATION OF COORDINATE 
SPACE 

Results of both numerical and analytical calculation 
presented in Fig. 3 reveal another unexpected feature of 
eigenfunction correlation which appears to be common 
to all ADRM. Surprisingly it is also present for the 3D 
Anderson model both in the metal and in the insula- 
tor phase (see Fig. 6). This is the negative eigenfunc- 
ton correlations for ui = \E — E'\ > Eq ^ b. In- 
deed, one can see from Fig. 3 and Fig. 6 that for large 
enough uj the correlation function C{uj) goes below the 
uncorrclatcd limit C{uj) — 1/N which corresponds to 

(|vl/,|2 Ivl/^f) = (|VI/,|2) (|xl,^.|2) ^ 1/^2^ dgj^o^g by E„ 

the value of uj where this limit is reached. For uj > Eq 
the correlation function C{uj) oc l/u;^ decreases down to 
zero. Such a behavior implies that two eigenfunctions 
separated by an energy difference to > Eq try to avoid 
each other. That is, if a site r is occupied in one of the 
states it should be predominantly empty in the other. 

To explain such a behavior the following cartoon is 
useful. Let us define a support of an eigenfunction as the 
manifold of sites {r} where |^'j(r)p is essentially non- 
zero. To construct such a support starting from a given 
site r with the on-site energy Sr we find all the sites in 
resonance with the site r, i.e. such sites r' which on- 
site energies Sr' obey the relationship \er' — Sr\ < \Hry\. 
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FIG. 5: Two-eigenfunction correlation function for the 3D 
Anderson model (orthogonal symmetry class) with a triangu- 
lar distribution of random on-site energies (solid symbols) and 
the critical PLBRM Eq.® with f5 = 1 and b = 0.42 (open 
symbols). The energy difference u = \E — E'\ is measured in 
units of mean level spacing. The insert shows the mean den- 
sity of states; the mobility edge corresponds to e = ±3.5. The 
energies E,E' were taken from the window (3.3,3.7) for the 
3D Anderson model and (-0.2,0.2) for the critical PLBRM. 
The slope of the critical power-law Eq.© is 0.52 in both cases 
which corresponds to d2/d = 0.48. 
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ing the conduction band states out of the on-site states in 
the tight-binding model. Then by construction the eigen- 
functions belonging to the same shell are well overlapping 
but those belonging to different shells do not overlap. 

From this cartoon it is clear that the physical meaning 
of the scale Eq is the width of the energy band corre- 
sponding to a single shell. Indeed, if the energy separa- 
tion Lu greatly exceeds the typical single shell bandwidth, 
the two eigenfunctions must belong to different shells and 
thus do not significantly overlap in space. On the con- 
trary, if LU is smaller than the single shell bandwidth, the 
two states typically belong to the same shell and thus 
overlap strongly no matter how sparse the shell support. 

The new energy scale Eq, which is the upper energy 
cut-off of the multifractal correlations, corresponds to a 
new length scale 
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FIG. 6: Eigenfunction avoiding for the PLBRM with 
b=0.42(stars), 3D Anderson insulator (red triangles) and 
metal (blue triangles). The dotted line corresponds to the 
limit of uncorrelated eigenfunctions; the solid line corresponds 
to the power law Points below the dotted line corre- 

spond to eigenfunction avoiding. 




FIG. 7: A cartoon of stratification of the coordinate space: 
different non-intersecting supports shown by different colours. 
Each support corresponds to a shell of states occupying this 
support and thus strongly overlapping; states belonging to 
different shells do not overlap. The stratification of space 
explains both strong correlations of states at energy separa- 
tion uj smaller that the single-shell bandwidth Eq and mutual 
avoiding of eigenstates for lo > Eq. 



Then the procedure should be repeated for all sites r' 
and so on. It is important that so obtained manifold {r} 
does not always include all the sites of the system. If this 
is the case, the whole coordinate space is stratified into a 
set of mutually non- intersecting supports (Fig. 7). 

Once the support is defined, one can build a shell of 
states on this support by making a linear combination of 
on-site states, pretty much in the same way as in build- 



which has a meaning of the minimum length scale of the 
fractal texture. In the c?-dimensional Anderson model the 
energy scale Eq can be estimated as 

E^r^V ^ D/Wc Dl{2d\n2d), 

where D is the total bandwidth. Estimating the DOS as 
p = l/{a^D) we find £o ~ aWl^'^^ where a is the lattice 
constant. 

Clearly the picture with a stratified coordinate space 
is possible for PLBRM Eq.([3]) with small enough h < 1 
when the single shell bandwidth Eq ^ b is small com- 
pared to the total bandwidth ^ 1. Amazingly, the 3D 
Anderson model which low-frequency critical features are 
well described by the critical PLBRM with b k, 0.42, also 
follows the predictions of the critical PLBRM for high 
frequencies co > Eq. This is a consequence of a rela- 
tively large value Wc = 16.5 of the critical disorder which 
results in Eq considerably smaller than the conduction 
bandwidth D. In particular its coordinate space must 
be stratified to explain the observed (see Fig. 6) mutual 
avoiding of eigenstates. 



VII. THE MULTIFRACTAL INSULATOR 

As has been demonstrated in Sec.V critical eigenfunc- 
tions with S, > L are strongly correlated in space. Here 
we consider the case of multifractal insulator where the 
localization radius ^ is large compared to relevant micro- 
scopic lengths (the lattice constant or elastic scattering 
length) but is much smaller than the system size L. We 
will identify a suitable random matrix model to describe 
this case and compare the properties of eigenfunction cor- 
relation in this model with those of the 3D Anderson 
model. 
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FIG. 8: Eigenfunction correlation in the ID Anderson insu- 
lator with rectangular distribution of on-site energies and pe- 
riodic boundary conditions. The disorder strength is Vl^ = 5 
(circles), Vl^ = 10 (stars). The inverse participation ratio is 
equal to 0.23 and 0.46, respectively. 



A. The ideal insulator limit 

We start by considering a limit of strong disorder when 
the localization length ^ ~ 1 and the multifractal nature 
of eigenstates does not show up. A common wisdom is 
that in the strongly localized regime the positions of the 
localization centers are completely uncorrelated. As it is 
shown in the Introduction, this leads to 

NC{lu) = 1, (19) 

which we will refer to as the ideal insulator limit. Fig.l 
shows how this limit is reached in the ensemble of banded 
random matrices. 

Note that C(a;)in this limit is much smaller than the 
self-overlap of |5'(r)p given by the inverse participation 
ratio Eq. ^ . Only for very small energy separations (typ- 
ically oc e"^/^) which we will not be considering here, the 
IPR limit can be approached. 

Now let us see how does the correlation function C{uj) 
look like for the strong Anderson insulator. The corre- 
sponding plot for ID Anderson model is shown in Fig. 8. 
It coincides almost exactly with the ideal insulator limit 
Eq.dm. 

The plot for the 3D Anderson model is shown in Fig. 9. 
On can see that NC{iS) is significantly enhanced at small 
energy separations and does not resemble at all the cor- 
relation function in ID Anderson insulator. 



FIG. 9: Eigenfunction correlation in the 3D Anderson insula- 
tor with rectangular distribution of on-site energies and peri- 
odic boundary conditions. The disorder strength is W = 80 
(green), = 60 (blue), W ^ AQ (red), W = 30 (purple). The 
system size is L = 20 for filled symbols and L = 8 for open 
symbols. The inverse participation ratio for the four insulat- 
ing systems is P2 = 0.72, 0.63, 0.44, 0.28 which corresponds 
to 4 = 1.0, 1.1, 1.2, 1.4 according to C = (9/47rP2)^/^ The 
change of the slope occurs at \E — E'\ — 5^. The slope for 
larger energy separations \E — E'\ > progressively increases 
with increasing W remaining smaller than 1. The insert shows 
the result for W — 60, L — 20 for the periodic (upper blue 
curve) and the hard wall (lower red curve) boundary condi- 
tions. 

B. Repulsion of centers of localization for (, 

In order to understand why the ideal insulator limit is 
not reached in the 3D case despite the ratio ^/L > 10 we 
compute numerically the probability distribution (PDF) 

F^{R) = {S{uj - E.^ + Era) 6{R - |r„ - r,„|)) (20) 

of the distance R = |r„ — rm| between the points r„ and 
r„i in real space {centers of localization) where |\['„^m(r)P 
has an absolute maximum, provided that the energy sep- 
aration between the states n, m is w. 

The results are shown in Fig. 10. It is seen that the 
function F^{R) is far from being independent of R (which 
would imply the lack of correlations between centers of 
localization). In fact, there is a repulsion of centers of 
localization at distances R < Rq ^ 12 which shows up in 
the decreasing probability density to find two centers of 
localization close to each other. Note that Rq is almost 
10 times larger than the localization radius ^ estimated 
from the inverse participation ratio P2- 
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be represented as follows: 




3D AM, W=60, L=20 



10=0.011 (black) 
10=0.0035 (blue) ■ 
(0=0.0015 (green) 



FIG. 10: The probability density Ea. (|20|l of having two cen- 
ters of localization at a distance R in real space and at a 
distance u in the energy space computed for the 3D Ander- 
son model in the strong localization regime {W = 60, ^ = 1.1, 
L — 20). The repulsive core i?o ~ 10 — 12 exceeds the "hard 
ball" limit 2^ by a factor of 5-6. 



An explanation to this fact of repulsion between cen- 
ters of localization rn and is based on the resonance 
interaction between states \l/„(r) and \l/m(r)if the energy 
distance between them is smaller than the typical overlap 
integral 14m (i?) oc e"^/^^. 

The size of the repulsion core i?o can be estimated from 
the equation: 



Vnm{Ro) — ^, 



Ro - 2^ In 



(21) 



The characteristic energy scale 6^ is the mean level sep- 
aration for states localized in the same volume .f'*. Thus 
the repulsion of centers of localization is a direct conse- 
quence of repulsion of energy levels for states confined 
in the same volume The energy scale depends on 
the strength of disorder and is of the order of the Fermi 
energy for strongly localized states. At ui the size 

of the repulsion core Rq may considerably exceed the lo- 
calization radius. 

The qualitative picture of repulsion of centers of local- 
ization can be quantitatively confirmed using the analyt- 
ical theory Eqs. pHll"4|) for the almost diagonal Gaussian 
RMT. To this end we look at the contribution of 3> 1 
to the sum in Eq. (|12p . Replacing the summation over n 
by integration we obtain the contribution to NC{lj): 



dn. 



(j(n)<w ^ 



(22) 



This equation can be easily interpreted using an elemen- 
tary perturbation theory. Indeed, for strongly localized 
states 5'm(rn) the eigenfunction correlation function can 



C„™ = ;^|*„(r)|2|^'„(r)|2 « (23) 

r 

« |*m(r„)P Wl' + l*«(rm)P E 

r r 

« |^'™(r„)P + |*„(r^)p. 

The amplitude at the tail of the wavefunction |^'n(rm)P 
with the maximum at a point r„ can be computed from 
the elementary perturbation theory in which the wave- 
function of the zero-th approximation corresponding to 
the energy £,„ is |^'i°^(r)p = 5r,r„: 



\^n{Vm)? = |«'m(r„)|' 



\Hnr. 



I) 



< 1. 



(24) 

The fluctuating on-site energy £„ is the main part of the 
eigenvalue En for a sufficiently strongly localized state. 
Thus we come to a conclusion that the amplitude of the 
wavefunction at a center of localization of the wave- 
function 5'„i is inversely proportional to (i?„ — i?„i)^ = 
and thus is strongly enhanced when uj <^ 5^. At the first 
glance this is in a contradiction with the common wis- 



dom that |^'„(r„)|" 



j|r„ r,„|/^ -^^ri-^jcj-^ jg apparently 



w-independent. The point is that the quantity |^'„(r)p 
has many accidental spikes due to resonances between 
on-site energies. The measure of such resonance points 
is small and for some (but not all) purposes one can ne- 
glect them and to approximate |\['„(r)p c>c el''"^'"™l/'>. 
The best known example when the two-spike eigenfunc- 
tion makes the main contribution is the low-frequency 
conductivity in the localized phase^^. As we will see be- 
low, here we deal with a very similar phenomenon. 

Now the correlation function C{uj) can be computed 
just by averaging over disorder and the distance R = 

I"?? I'm . 



C{uj) = I FUR){Cnm{R))d''R (25) 



FUR)' 



where F^{R) is the PDF defined by Eq.pD|. 

Comparing Eq. ((25l) for d = 1 with Eq. (f22|) we see that: 



FUR) = l/N, (i?»i?o), 



(26) 



where Rq is found from the condition (t{Rq) = lo similar 
to Eq. (|2T|) . 

In the opposite limit i? <C i?o 7 or w <C (t{R), we have a 
resonance enhancement [^^(rn)^ « 1/2 and C„m ~ 1- 
Then the comparison of Eg. ([^5]) with Eq. p^ yields: 



FUR) 



3iV \a{R) 



(i?«i?o). (27) 
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can be formally taken into account in the random ma- 
trix formalism Eg. ifT^ if one assumes the follo"wing re- 
lationship bet"ween the d-dimensional vector R and the 
difference of matrix indices n ~ m: 



d{n- m) ^ VLdR'^-^ dR, \n - m 



^i?^ (29) 



■where fid is the total solid angle in the d-dimensional 
space. 

In particular Eq.f^H) suggests that for exponential lo- 
calization the correct truncating factor in Eq. ([5]) has the 
form: 



exp 



I I \ ^/d-' 

\n — to| \ 



B 



This sets the exponent 77 in Eq.([5]) equal to: 



(30) 



(31) 



FIG. 11: A comparison of eigenfunction correlation functions 
for 3D Anderson model with W = 30, L = 20 and the trun- 
cated critical RMT Eq.(7) with b = 0.42, B = 5, ?? = 1/3 
and iV = 2000. The IPR takes values of P2 = 0.28 and 0.25, 
respectively. The scale of tj = \E — E'\ is different in those 
two cases by approximately a factor of 11. 



C. Logarithmic enhancement of correlations of 
locahzed eigenfunctions and the Truncated Critical 
RMT 



No-w let us consider Eq. ([25)l for d > 1 assuming that 
all states are exponentially localized and thus (|i?nmP) oc 
exp[— i?/^]. We also assume for simpHcity that F^{R) = 
N-^ 9{R ~ Ro), where Rq = 2(,\n{S^/ujj. Then one im- 
mediately obtains from Eg. ips)) that due to the phase 
volume factor R'^^^ the correlation of exponentially lo- 
calized eigenfunctions depends crucially on the dimen- 
sionality of space. Namely, for d = 1 the ideal insulator 
limit Eg. (fin)) is reached for sufficiently small ^ (see Fig.l), 
"while for d > 1 and uj <^ the correlation function ac- 
quires a logarithmic in uj enhancement factor: 



NC{uj) - C 



d-d. 



(28) 



The physics behind this result is similar to the one 
■which leads to the selebrated Mott's la^w^iS^ (^{^) oc 
ln'''~^^{S^/uj) for the ac conductivity a{uj) at w ^ S^. 
The difference is that the contribution to conductivity 
from the resonance states -with the distance R between 
the points of maximal amplitude is proportional to the 
square of the dipole moment d"^ (x R^, so that the phase 
volume factor i?'^"^ gets multiplied by R^ resulting in 
emergence of the logarithmic factor ln'^^^(5^/(jj) instead 
of ln'^~^((5^/(jj) in our case. 

Below we obtain this result for the truncated critical 
RMT defined by Eq.dg. The phase volume factor R'^'^ 



Then Eq. (jl2p can be used which is convenient to rewrite 
in the following form: 



NC{uj) « k{uj) ^ - 



Ji^y) 

3;j(lnCT2(n))|„^„(y) 



(32) 



where n{y) is found from the equation 

,2 



y 



UJ 



and 



fiy) = 



Vv 



dt— = { ! 
y 



772, 2;>i 

(33) 

For the truncated critical RMT Eq. ([8]) with rj = 1/d one 
finds: 

1 n 



^(lna2(n)) 2+^^)'^'' 



where 



n{y) 



(34) 



The integral in Eq. (|32|) is well convergent and thus 
mainly contributed by y ~ 1. This makes it possible to 
obtain a simple analytical expression for NC{uj): 



NC{uj) 
where 



'cdS In'^-i (^) , u:<^5^^b/B 



Co 



(35) 



Cd = 2 



d I f{y)dy^2''d. 
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The first line of Eq.(l35l)is consistent witli Eq.lUHl) in 
which d2/d ~ 6 < 1 and ^ - B^/<^. The power-law 
behavior in the second line of Eg. ([35]) is a remnant of 
the critical behavior Eq.®. It exists only for consider- 
ably large B ^ 1 where (5^ ^ i?o = cofe, i.e. only in the 
multifractal insulator. For a very strong insulator with 
6^ 3> -Eq J the localization radius is smaller than the min- 
imal length scale £o of the fractal texture. This is the 
region of an ordinary insulator where the entire localiza- 
tion volume is more or less homogeneously filled. 

Thus the eigenfunction correlation function for the 
truncated critical RMT Eq.([8])interpolates between the 
behavior given by Eq. ([28|) at a; <C (5^ (or = 
(pcj)^^/'* ^ ^) and the critical behavior Eq.® which 
is valid for £o <^ L^j <^ At ~ ^ both asymptotic 
forms are apparently matching with each other. 

The physical picture that leads to such a behavior is 
the following. There are two distinct regions < ^ 
and > ^ where physics of eigenfunction correlations 
is entirely different. In the case < the characteris- 
tic length L^^ has a meaning of the period of beating in 
the overlap of two fractal eigenfunctions inside the local- 
ization volume. The regions where two fractal supports 
match well with each other alternate with the regions 
with a strong mismatch between them, very much like in 
the case of two grids with slightly different periods. The 
regions of strong overlap make the main contribution 

to the eigenfunction correlation function C{ijj) which is 
of the order of the IPR of a multifractal metal (see next 
Section for more details) with the system size equal to 
the localization radius ^ and the correlation length equal 
to the size of the well overlapping regions. To ob- 
tain the correlation function C(w) one has to multiply 
Cm(w) by the probability for the entire localization vol- 
umes to overlap. This probability is £f /L'^, as for oj > 
there is no correlations in the positions of the localization 
volumes. Thus we obtain the critical power-law Eq.®: 
NC{uj) ^ {L^ll^f-''^ = {E^lu^f-'^^l''. 

For > ^ physics of eigenfunction correlations 
changes drastically. Now localization volumes are sta- 
tistically repelling each other and the overlap is only due 
to the tails. In this region the length scale loses its 
physical meaning which is taken over by the length scale 
i?o given by Eq.([2T|). 

The overall shape of C(lo) with the logarithmic en- 
hancement factor ln'^^^((5^/ti;) obtained within the Trun- 
cated Critical RMT describes the numerical results on 
the 3D Anderson model very well (see Fig. 11). The ab- 
sence of this factor at d = 1 explains the qualitative dif- 
ference between the case d=\ (see Fig. 8) and d — Z (see 
Fig. 9). This difference is essentially due to a competition 
between two effects (i) repulsion of centers of localization 
and (ii) resonance enhancement of overlap by tails. The 



first effect tends to decrease the probability of the overlap 
of localization volumes. The second effect increases the 
eigenfunction overlap by means of tails. In the ID case 
these two effects compensate each other and the result 
is the same as one would obtain for completely uncor- 
related positions of localization volumes and the typical 
exponentially decreasing tails. In higher dimensions the 
enhancement of overlap in the tail region prevails because 
of the increased volume of those regions. 

Concluding this subsection we claim that the truncated 
critical RMT provides an excellent description of the 3D 
Anderson insulator both in the strong localization re- 
gion (see Fig. 9) and in the region of multifractal insulator 
where the localization radius is large and the correspond- 
ing scale 5^ is small compared with the upper cutoff Eq 
of multifractal correlations. Because of the limited size 
of the 3D lattice this latter region is out of reach for nu- 
merical simulations on the 3D Anderson model, and the 
random matrix theory is the only mathematical model 
which properly describes physics of the multifractal in- 
sulator. 



D. Super-critical PLBRM 

Note that there is another RMT Eq.® suggested in 
Ref.™ as a candidate to describe eigenfunction correla- 
tions in the multifractal insulator. Below we show that 
this super- critical PLBRM is principally flawed, as it cor- 
responds to a power-law localization which is not the case 
in the 3D Anderson model. 

This can be best demonstrated by Ea. ((24|) in which 
\Hnm\ ^ cr(-R) = (b/R)" with a > 1. Accordingly, the 
typical scale for the repulsion of centers of localization is: 



In Fig. 12 we plot the results of numerical calculation of 
the PDF F^{R) for the super-critical PLBRM Eq.®. 
The characteristic scale R^j where F^{R) reaches its max- 
imum, is well seen in this plot. 

The analytical treatment based on Eg. ([5^ yields for 
this model: 

j_ 

NC{lo) - (^^^ " ^ R^, (37) 

where Ea — (cq6)", and 

7r3/2 1 

As well as the entire approach based on Eq. (fT2|) . the 
above results are valid when b is the smallest relevant 
parameter. In the problem of PLBRM with a close to 
1, there is a competition between the small parameters 
|1 — a| and b, so that the validity of Eqs. ((5^ . ([57|) requires 
also a — 1 ^ b. 
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FIG. 12: The correlation function of centers of localization 
F^{R) for the super-critical PLBRM Eq.® with a = 1.1, 
b = 0.1, iV = 1000 for Lu = 0.011 (black), 0.0035 (blue), 
0.0015 (green). The insert shows the ti;-dependence of Rui 
where Fuj{R) reaches its maximum. The finite slope of Fui{R) 
for R ^ Ru, is a finite size effect which was neglected in 
Eq.JleJ. 
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FIG. 13: Two-eigenfunction correlation for the super-critical 
PLBRM Eq.(|9l) with a — 1 >> b calculated analytically using 
Eas. (llH14|l . The power law C{uj) oc tj"^/" is vahd for all 
energy separations corresponding to 1 < R^^ < N. The onset 
of the plateau moves to a; —> in the limit N ^ ex. The 
ideal insulator limit is reached by decreasing the slope with 
increasing a. 



We see that in the infinite system N ^ oo the power- 
law Eq. (j37p in C(a;) is not restricted at small lu, and no 
energy scale similar to 6^ emerges. This can be explained 
only if we assume that the localization length for a—1 >> 
b is of order one. Then for all energy separations u ^ Ea 
the repulsion core Ruj ^ and no qualitative change in 
the correlation function occurs until hits the system 
size N . For smaller uj the correlation function is almost 
a constant. This quantitative analysis is illustrated by 
Fig.13. 

The region of a that could describe the multifractal 
insulator with large ^ 1 corresponds to a — 1 ^ 6. In 
this case an energy scale similar to (5^ should appear. It 
can be found from the condition 

RM=5,=aa), ^6^^^. (38) 

At \E — E'\ = (5^ the slope on the log-log plot of C{lj) 
should change from the critical value at (5^ < |£' — < 
Eq to a different (but constant) a-dependent value at 
\E ~ E'\ < S^. This change of the slope is clearly seen in 
the numerical simulations on the super-critical PLBRM 
presented in Fig. 14. It appears that in all cases studied 
the slope ^ at \E ~ E'\ < is larger than that at |i? — 
E'\ > 5^. This is in a clear contradiction with the results 
(see Fig. 9) obtained in the 3D Anderson insulator. 

An important conclusion we can draw from the above 
analysis of the super-critical PLBRM is that the corre- 
lation function C(u;) for uj < is the power-law in this 
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FIG. 14: Two-eigenfunction correlation for the super- critical 
PLBRM Eq.®. Numerical results for a 1 < 6 (this plot 
corresponds to a = 1.15, b = 0.45, A'^ = 5000) show that 
the exponent /i of the power-law C(tj) cx changes at 

\E — E'\ = with the larger value corresponding to smaller 
energy separations. 



model. This can be traced back to the power-law char- 
acter of localization in the super-critical PLBRM which 
is not the case in the disordered lattice models (such as 
the 3D Anderson model) with short-range hopping inte- 
grals. This is the reason why the super-critical PLBRM 
is not suitable to describe the insulating phase of the 3D 
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Anderson model. 



VIII. SEARCH FOR RANDOM-MATRIX 
MODEL FOR A MULTIFRACTAL METAL 

A. Anti-truncated critical RMT 

Surprisingly, the natural counterpart to the truncated 
critical RMT Eq.® which is defined by Eq.dTO]) ("anti- 
truncated" critical RMT) does not describe extended 
states in the multifractal metal. The reason is that this 
model possesses two low-frequency system-size indepen- 
dent energy scales instead of the single scale which is 
associated with the size ^ of a multifractal cell in Fig. 2c. 
In order to see this we analyze the analytical formulae 
Eas. pim4p with the variance defined by Eq. (jlOp . 

To this end we expand the summand of Eq.((T2]) in 
l/w <C 1 to arrive at the formula similar to Eq. (|^^ but 
with the upper limit of integration equal to the correla- 
tion radius ^: 



2 

NC[i.) = - / 



(T^(n) dn. 



(39) 



where ^ = B/2h according to Eg. p^ . 

Substituting Eq. fTO)) for a{n) we arrive at: 



NC{uj) 



62 




■ 1 


2 " 


Cj2 


f 

Jb/uj 







h 

UJ 



2^ 



tj2B2 



(40) 



Eq. (j40|) is valid at lo > b/^ when the upper limit of in- 
tegration is larger than the lower limit. This sets the 
energy scale wi = b/£, = 2b^/B. 

Another scale ^2 = 1/ (26^) gives the cross-over scale 
that separates the critical 1/uj behavior and the l/w^ 
behavior that takes place for indermediate frequencies 
uji < oj < UJ2- While the scale cui — b/^ (similar to 
the scale 6^ = b/B in a ID insulator) determines the 
onset of the low-frequency plateau, the second relevant 
scale UJ2 that appears in the model Eq. PU)) seems to have 
no physical meaning. Indeed, the existence of this scale 
leads to a characteristic form of the correlator C{oj) which 
log-log plot has a significant slope increase just before it 
drops to zero at the plateau (see Figs. 15, 16). 

We did not find the behavior of such type in the 3D 
Anderson metal (see Fig. 17). The plot in Fig. 17 clearly 
shows a saturation' 
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Cito) 



at < ti>i: 



tP2 = A„ 



Si" 



0.5. 



(41) 



However, there is no evidence of a maximum in the slope 
just above the onset of the plateau. 



B. Sub-critical PLRBM 

Now we consider the sub-critical PLBRM ensemble de- 
fined by Eq.® with l/2<a<l. In this case analytical 
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FIG. 15: Numerics on the anti-truncated critical RMT of the 
orthogonal symmetry class. The predicted analytically non- 
monotonous behavior of log-log slope controlled by the energy 
scale u!2 is well seen both for small and large system sizes. The 
solid line is a fit according to Eq. (|40)) . 




FIG. 16: Eigenfunction correlation for the anti-truncated crit- 
ical RMT Eq.dlJ (red curve) and for the sub-critical PLBRM 
Eq.® (blue curve) of unitary symmetry class computed an- 
alytically from Eas. (|llj | T4|l . For the anti-truncated critical 
RMT the non- monotonous behavior of log-log slope is similar 
to the one in Fig. 15. 



arguments similar to Eqs. (|39l40p predict only one rel- 
evant energy scale uji = 6/^" such that for oj < uii the 
correlation function C{uj) is constant and for lu > lui (but 
Lo < Eq ^ b°') it is a, pure power law C{uj) oc lu^~ . Thus 
the sub-critical PLBRM is free from the drawback related 
with the unphysical second energy scale. For comparison 
we plotted the analytical results for the anti-truncated 
critical RMT and for the sub-critical PLBRM in Fig. 16. 
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FIG. 17: Eigenfunction correlation in the 3D Anderson model: 
extended states. The disorder strength is W = 2 (green), 
W = 5 (red), W = 10 (blue). The system size is L = 8 (di- 
amonds) L = 16 (circles) and L — 20 (squares). The corre- 
lation length is estimated from the inverse participation ratio 
as follows C = (P2iVMm)^/<'*"'*'' (^m = 0.5, (^2 = 1.3) and is 
equal to: C = 5.8, 3.2, 2.8 at W = 10, 5, 2, respectively. In 
the insert we compare results for W — 5, L — 20 for the peri- 
odic (red filled squares) and the hard-wall (red open squares) 
boundary conditions. For hard-wall boundary conditions the 
correlation function looks "more critical", as the critical point 
Wc = 15.2 in this case is closer to W = 5. 




FIG. 18: Eigenfunction correlations in the sub-critical 
PLBRM of the orthogonal symmetry class. Left panel (a): the 
small-6 limit b = 0.001, a = 0.8 at A?" = 1000 (squares) and 
N — 200 (diamonds). The critical correlations for b = 0.001, 
a — 1 and A'^ = 1000 are shown by stars. The solid line cor- 
responds to the power-law u)~^. The slope of the sub-critical 
curve is larger than 1 (correlations are short-range in the en- 
ergy space). This case is relevant for the metallic phase of the 
Anderson model in very high dimensions. Right panel (b): the 
case b = 0.4, a = 0.95 at A = 1000 (squares) and A = 200 
(diamonds) is relevant for the Anderson model in d = 3 ("mul- 
tifractal" metal). The slope is less than the critical (which in 
turn is less than 1) and is almost size-independent. Correla- 
tions are long-range in the energy space. 



It is seen that the overall shape of the blue curve for 
sub-critical PLBRM is much closer to the results of 3D 
Anderson model of Fig. 17. 

Note however, that the power-law emerging in the an- 
alytical results for the sub-critical PLBRM has an ex- 
ponent 1/a which is larger than the critical exponent. 
Computer simulations (see Fig. 18a) on the sub-critical 
PLBRM with very small b confirm this analytical result 
as cxD extrapolation and show that the slope in- 

creases with increasing the system size N. However, the 
slope of the corresponding curves for the 3D Anderson 
model of Fig. 17 is almost independent of the system size 
and is equal or smaller than the critical slope. 

The reason for the discrepancy is that the analytical 
result for the slope ^ = 1/a corresponds to the limit 
6^0. At a finite b the slope decreases with increasing 
b and at a sufficiently large b may become smaller than 
the critical one. It is reasonable to assume that at small 
1 — a and b this happens at 6 ~ (1 — a). The relevance of 
the parameter b/{l — a) is also seen from the expression 



for the correlation length ^ 
^ - exp 



1 - a 



, aa = In { ) + const. (42) 



which was found (up to a constant of order one) analyti- 
cally from Eq. (|16p . This expression is apparently mean- 
ingless for 6" 3> (1 — a) where Ca may become negative. 

Numerical simulations on the sub-critical PLBRM 
with 1 ~ a < b (e.g. for a — 0.95 and b — 0.4 rel- 
evant for the 3D Anderson model) show (see Fig. 18b) 
that the log-log slope of C{uj) is somewhat smaller than 
the critical one and is almost independent of the matrix 
size N. Thus the sub-critical PLBRM shows exactly the 
same character of eigenfunction correlations as in the 3D 
Anderson metal (see Fig. 17). 

Two parameters of the sub-critical PLBRM allow to 
simulate the effect of the finite correlation length (choice 
of a < 1) and the dimensionality of space (choice of b). 
Note in this connection that for the disorder strength W 
significantly smaller than the critical value Wc, not only 
a but also b is M^-dependent. The point is that in the 
3D Anderson model the variance of the on-site energies 
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fluctuations is proportional to W"^, while the off-diagonal 
hopping integral is equal to 1. This implies that the ra- 
tio of a typical off-diagonal to a typical diagonal elements 
controlled in Eq.® by the parameter b should scale like 
1/W. As the log- log slope of C{uj) decreases with in- 
creasing b, moving away from the Anderson transition 
into the metallic phase W < Wc has an effect of decreas- 
ing the slope. On the insulator side of the transition the 
situation is opposite and one should expect an increase 
of the slope (for 5^ < uj < Eq) with increasing W. Fig. 9 
shows that it is apparently the case. 

Another relevant note is that for the Anderson model 
in higher dimensions the correlation dimension d2 de- 
creases. This can be modeled by a decreasing parame- 
ter b. Then the analogy with the sub-critical PLBRM 
suggests that for sufficiently high dimensions d > dc the 
behavior in the d-dimensional Anderson model should be- 
come similar to the one in Fig. 18a. Namely, the exponent 
fi in Eq.([7]) may become larger than 1. This changes qual- 
itatively the eigenfunction correlations, as they become 
effectively short-range in the energy space. In particu- 
lar, the return probability" which is proportional to the 
Fourier transform of C{ljj) behaves in the time interval 
K/Eq < t < h/5(^ as P{t) cx ^-(1-'') for yit < 1 and is a 
constant for /i > 1. 

We believe that this qualitative change in the eigen- 
function statistics (if confirmed for a d-dimensional An- 
derson model with d > dc) should lead to dramatic phys- 
ical consequences marking a transition to a new metallic 
state. 



IX. CONCLUSION 

In conclusion we list the main results obtained above. 
The most important of them is the persistence - beyond 
the point of localization transition- of the critical power- 
law in the dependence of the eigenfunction correlation 
function C{uj) on the energy separation lo and the re- 



lated enhancement of C{lo) at (5^ <C o; <C i?0: where 
is the mean level spacing in the localization/correlation 
volume and Eq is the upper energy cut-off of multifrac- 
tality. This enhancement leads to an enhancement of 
matrix elements of local electron interaction which may 
result in, e.g. an enhancement of the superconducting 
transition temperature in the vicinity of the Anderson 
localization transition^^. Another important observation 
is that the enhancement of correlations a.i lo < Eq is al- 
ways accompanied by the depression at cj > _Eo, both 
phenomena being the consequences of the stratification 
of the coordinate space into mutually avoiding supports 
of the fractal structure with well overlapping eigenfunc- 
tions living on each of them. An independent - but 
also important- phenomenon is the logarithmic enhance- 
ment of C((jj) in the 2D and 3D Anderson insulator at 
w < (5^ (and the absence of such enhancement in the 
quasi- ID disordered wire). It is a result of a competi- 
tion of two simultaneous phenomena: the repulsion of 
centers of localization and the resonance enhancement 
of the eigenfunction overlap by tails. Both phenomena 
are studied quantitatively within the Truncated Critical 
Random Matrix model which is suggested as a univer- 
sal tool to describe the localized eigenfunctions with a 
multifractal texture. We also show that the sub-critical 
Power-Law Banded Random Matrix Ensemble suggested 
in Rcf.^"* describes the multifractal metal reasonably well. 
From the analytical solution for this RMT we conclude 
that a critical dimensionality dc may exist above which 
the d-dimensional Anderson model has an unusual metal 
phase characterized by an effectively short-range correla- 
tion function C{ljj) cx lo~^ with > 1. 
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